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ABSTRACT 

We focus on Hinode Solar Optical Telescope (SOT) calcium II H-line observations of a solar quiescent 
prominence (QP) that exhibits highly variable dynamics suggestive of turbulence. These images capture a 
sufficient range of scales spatially (^^0. 1-100 arc seconds) and temporally (~16.8 s - 4.5 hrs) to allow the 
application of statistical methods used to quantify finite range fluid turbulence. We present the first such 
application of these techniques to the spatial intensity field of a long lived solar prominence. Fully evolved 
inertial range turbulence in an infinite medium exhibits multifractal scale invariance in the statistics of its 
fluctuations, seen as power law power spectra and as scaling of the higher order moments (structure functions) 
of fluctuations which have non-Gaussian statistics; fluctuations Sl(r,L) = I{r+L)-I(r) on length scale L along 
a given direction in observed spatial field / have moments that scale as < 6l{r,Ly' >r^ L'^^p\ For turbulence in 
a system that is of finite size, or that is not fully developed, one anticipates a generalized scale invariance or 
extended self-similarity (ESS) < 5l{r,Ly' >^ G(L)^*'^'. For these QP intensity measurements we find scaling 
in the power spectra and ESS. We find that the fluctuation statistics are non-Gaussian and we use ESS to obtain 
ratios of the scaling exponents ({p): these are consistent with a multifractal field and show distinct values for 
directions longitudinal and transverse to the bulk (driving) flow. Thus, the intensity fluctuations of the QP 
exhibit statistical properties consistent with an underlying turbulent flow. 

Subject headings: Sun: corona - Sun: prominences - magnetohydrodynamics (MHD) - plasmas - turbulence 



1. INTRODUCTION 

Solar prominences or filaments in the lower solar corona 
are relatively cool, dense plasma structures with temperatures 
of about 10"* K. Solar filaments can be seen on the disk, 
whilst prominences are observed above the solar limb. In 
practice, they are classified in three main categories accord- 
ing to their location on the Sun, namely active, intermediate 
and quiescent. The latter usually occur on the quiet Sun at 
high latitudes and as a consequence are also known as "polar 
crown" prominences, while active and intermediate filaments 
are oft en observed at low latitudes associated with active re- 
gions (lEngvoldlll998!) . All prominences originate from fila- 
ment channels and develop above the polarity inversion line. 
They show many different morphologies and dynamics (see 
iMackav et al. 2010, for a recent review). 

The Hinode Solar Optical Telescope (SOT) provides ob- 
servations of solar prominences revealing detailed internal 
dynamics at unprecedented spatio-temporal resolution. In 
particular, dynamics associated with quiescent prominences 
(QPs) are seen to exhibi t spatio-temporal ev olution character- 
ized by high variabihty dBerger et alj|2008l) . Most of the QPs 
in the SOT dataset appear vertically structured and dominated 
by upward and downward transport of matter; the ascend- 
ing flows appear dark and are faster (25 kms~'on average) 
than the descending flows (about 10 kms~'). The upflows 
have often been obser ved to ultimately evolve into vortices 
jLiggett & ZirinI [19841) and are considered to be associated 
with small scale turbulence dBerger et al.ll20ldl) . The promi- 
nence is a low-/? plasma with electron density ~ 10^'cm"^ 
and temperatures up to ^ 10"* K dTandberg-Hansseiil 1 1 995h ; 
these typical parameters suggest that the upward flow is su- 



personi c. Indeed, evidence of bow-shock compressions are 
seen in iBerger et alj d2010l) and the corresponding Reynolds 
number is estimated as ~ 10^. The question then immediately 
arises as to whether the observed fluctuations do in fact corre- 
spond to a turbulent flow. 

Many models have been developed to describe possi- 
ble scenarios for the production of dynamical structures 
in the corona. The local magnetic field is suggested 
to play a key role as it is thought to b e the driver of 
the prominence th r eads (e.g.. Low & Hun dhausenI 119951 ; 
iFoullon et ai] 120091; iHershaw' et all |2011). Recently, in 
strongly inhomogeneous coronal plasma structures, processes 
such as magneto-t hermal convection in solar prominences 
dBerger et al."201 V) and Kelvin-Helmholtz instabilities in the 
corona (Foullon et al. 2011) have been suggested as mech- 
anisms for the generation of dynamical structures. Fur- 
thermore, the observed complexity of the coronal mag- 
netic field may be generated by photospheric turbulence 
dAbramenko et al.ll2008t iDimitropoulou et alJl20"o'9) . Intrigu- 
ingly, correlations between outer corona and solar wind have 
also been found in the statistics of large-scale density fluc- 
tuations (Telloni et al. 2009) suggestive that the signature of 
coronal turbulence is convecte d with the solar wind plasma 
dMatthaeus & GoldsteinI 119861) . To distinguish these pro- 
cesses from turbulence evolving in-situ (locally) in the flow, 
we will apply analysis methods that have been specifically de- 
veloped to quantify finite range fluid turbulence. 

The characteristic, reproducible properties of a turbulent 
flow are statistical in nature. They characterize a scale in- 
variance of the statistical properties of fluctuations - that is, 
these properties are unchanged as we move from scale to 
scale subject to a rescaling. Thus in a fully evolved magneto- 
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hydrodynamic (MHD) turbulent flow in an infinite medium, 
one finds power law dependence of the physical observables 
of the flow - e.g., the velocity and magnetic field fluctua- 
tions have power law power spectra over a range of scales, 
which is identified as the inertial range of the turbulence. 
As po wer law power spectra are not unique to turbulence 
jSorne tte 2000) an d do not uniquely cha racterize the scaling 
of the fluctuations dChapman et al.ll2008l) . multifractal scaling 
of the higher order moments (structure functions) of fluctua - 
tions is also needed to identify a turbulent flow (Frisctf 19951) . 
These statistical methods have been applied extensively to in- 
situ observations of the outgoing flow from the solar corona, 
namely the solar wind, and have established its turbulent 
character (e.g., [H orbury & Balogh 1997; Sorriso-Val vo et alj 
[19991: iPagel & Balogh,2003>: ,Bruno & C arbone 2005i). How- 
ever, these in-situ observations are typically singl e or from a 
few point in space so that Taylor's hypothesis (Tayloi^ ll938l) 
is usually evoked to characterize the scaling properties of the 
flow. Here, SOT observations of a QP provide a direct obser- 
vation of the spatial field of fluctuations. 

In this paper we present the first application of these sta- 
tistical methods to the spatial intensity field of a long lived 
solar prominence. The SOT images are of intensity, rather 
than velocity, and so intrinsic to our analysis is the assump- 
tion that the moving structures in the images follow the 
flow, acting as markers or passive scalars for the plasma 
dynamics (e.g. the intensity measurements can be treated 
as proportional to the squared density of the plasma flow). 
This assumption is supported by the correspondence be- 
tween traceable UV motions and true mass motions, found 
using combined im aging and Doppler data of prominences 
jKucera et alj|2003h . Importantly, the QP is of finite physical 
size, and the turbulence may not be fully developed. Under 
these circumstances, we anticipate a generalized form of scale 
invariance, that is, generalized similarity (scale invariance) 
also known as Extended Self-Similarity (ESS) (iBenzi et al.l 
11993b . Generalized simila rity has been seen in the fast solar 
wind (C arbone et al. 1996: lHnat et aI.l2005l:lNlc"ol et alj2008l: 
IChapman & Nicol 2009), in laboratory simulations of MHD 
turbulence (Dudsonetal. 2005; Dendv & Chapman 2006) 
and i n hydrodynamics (Grossmann et al.l ll994t IBershadskiil 
l2007h . We find that the intensity fluctuations in the QP do in- 
deed exhibit quantitative features consistent with a finite size 
turbulent flow, namely, ESS, multifractality and non-Gaussian 
statistics. 

2. THEDATASET 

The Hinode spacecraft was launched in September 2006 
and moves in a sun-synchronous orbit over the day/night ter- 
minator, allowing near-continuous observations of the Sun. 
The SOT on board Hinode is a diffraction-limited Gregorian 
telescope with a 0.5 m aperture, which is able to provide im- 
ages of the Sun with an unprecedented resolution up to 0.2 
arcsec and cadences between 15 and 30 sec. The Broadband 
Filter Imager (BFI), one of the four instruments of the Fo- 
cal Plane Package on the SOT, provides observations over 
a range of wavelengths (380-670 nm) which distinguish dif- 
ferent coronal structures. We use the Ca II H spectral line 
(396.85 nm) images of a QP observed by the SOT on the 
north-west solar limb (90W 52N) on November 30th, 2006 
(s ee the time e volution of the QP of interest in animation 1 
of iBerger et al.l (2008)). The time interval considered covers 
-4.5 hrs, from 01:00:00 UT to 05:30:00 UT corresponding 
to about 1000 images with a cadence Af = 16.8 sec on av- 



erage at a spatial resolution of 0.10896 arcsec per pixel, that 
is, one pixel corresponds to Ar — 77.22 km on the solar sur- 
face; each image is 800 x 420 pixels. The images have been 
calibrated (normalized to the exposure time) and aligned with 
respect to the solar limb. Furthermore, these specific observa- 
tions are along a line of sight that is to a good approximation 
perpendicular to the prominence sheet. 

Figure [T] shows the first frame of the dataset. Note the dif- 
ferent structures: large scale structures appear brighter at the 
edge of the prominence while at smaller scales, bright and 
dark threads alternate within the plasma sheet. We will ex- 
amine fluctuations in space by taking differences in intensity 
along directions longitudinal (vertical) and transverse (hori- 
zontal) to the direction of upward/downward flow. This proce- 
dure is shown by the overlaid grid which is made of 10 strips 
labelled as strips LI to L5 along the longitudinal direction and 
strips Tl to T5 along the transverse direction; each strip is 10 
pixels wide. We will also examine fluctuations in time, that 
is, from one image to the next. Five white squares, labelled A 
to E, with size 21x21 pixels, indicate the regions over which 
the respective intensity time series are formed across all the 
images. 

In order to improve statistics we will construct local spa- 
tial a verages and will pres ent the variation about these aver- 
ages (iDudok de Witil2004h . The procedure used to analyse 
the intensity measurements in the strips consists in calculat- 
ing statistical quantities for small ensembles of 10 neighbour- 
ing rows (columns) for each strip along the horizontal (verti- 
cal) direction and then performing an average across the strip 
width. For example, the mean value of the intensities for strip 
Tl will be the average over the mean values calculated for 
each of the 10 rows within Tl . The same procedure is adopted 
for the analysis in the time domain: the statistical quantities 
calculated for each time series associated with the pixels that 
compose the square are averaged over the 21x21 pixels. 

Figure |2] plots the variation in intensity I(r) for strip T5 
versus pixel position r (top-left panel) and /(f) for square D 
versus time (bottom-left panel). They fluctuate strongly in 
their first differences, which are defined in space as 6l{r, L) = 
I(r+L)-I(r) with L = 1 pixel (top-right panel) and in time as 
(5/(f,T) = I(t + T)-I{t) with T = At = 16.8 sec (bottom-right 
panel). This is a typical aspect of a stoch astic process includ- 
ing turbulence jKantz & Schreiberil 19971) . 

3. POWER SPECTRA 

Fully evolved MHD turbulence is self-similar in character 
exhibiting power law scaling in the power spectrum. We thus 
analyse the power spectral densities (PSDs) of the intensity 
measurements both in the time and space domains. 

The left panel in Figure [3] shows the PSDs of the intensity 
measurements for strips T4 and T5 along the transverse direc- 
tion and LI and L2 along the longitudinal one. All power 
spectra are dominated by two main slopes: at small wave 
numbers the spectra scale as ~ consistent with a B rown- 
ian process, that is, additive noise (iPercival & Waldenll200Q) . 
while at larger wave numbers the spectra scale as ~ kr°' with 
spectral index a reported in TableQjfor each strip and sugges- 
tive of non-trivial dynamics. The a values are estimated by 
extracting the gradient of the linear fits to the plots in Figure 
[3] (left) within the wave number ranges 2.43-4.55 Mm"' for 
the longitudinal strips and 2.07-4.61 Mm"' for the transverse 
strips. The a values found are distinct from -5/3, which is 
the value expected for the Kolmogorov spectrum for an ideal 
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turbulent flow. This is not surprising since these observations 
are integrated, line of sight intensity measurements. However, 
we should still expect line of sight measurements to capture 
qualitative features of turbulence (such as non-Gaussian fluc- 
tuations, multifractal scaling and ESS) whilst not necessary 
giving the same numerical values of scaling exponents as in- 
situ point observations. 

The time series associated with squares A to E reveal dif- 
ferent dynamics: they have power law power spectra in the 
frequency domain with fitted spectral indices a in the fre- 
quency range 1-20 mHz very close to -1 (see Table [T] and 
the right panel of Figure [3]l. This ^ 1 // scaling may be 
simply attributable to a "random telegraph" process, that 
is, how a series of uncorrected pulses or features in the 
flow movin g through the line of sight of the observations 
(iKaulakys & Mes Kauskas 1998; Kaulakys et al. 2005). We 
estimate the "maximum observable speed" of structures mov- 
ing past a line of sight as m = Ar/Af ^ 4.6 kms~'. Since the 
prominence flow has a bulk velocity (m fio^ ^ 25 km s~' ) larger 
than u then, at a given pixel, intensity fluctuations are moving 
too fast for us to observe correlations in time. In other words, 
the time needed to catch a coherent structure (e.g., up-flows), 
at fixed space coordinates across two consecutive frames, is 
much shorter than the cadence, therefore, all the moving flows 
in the prominence appear decorrelated in time. 

4. PROBABILITY DISTRIBUTION 

We now investigate the statistics of the intensity fluctua- 
tions in the space domain, 5l{r,L) = I(r+L)-I(r) with length 
scale L, and in the time domain, 6l(t,T) = I(t + T)-I{t), with 
time scale r. Turbulent fluctuations in the inertial range in- 
variably possess a non-Gaussian "heavy tailed" probability 
density function (PDF) that arises from the intermittent na- 
ture of the energy cascade in the flow (iMarsch & Tul Il997t 
iSorriso-Valvo et al.l[T999HHnat et al.ll200 2'K 

The left panel in Figure |4] shows the PDF of the intensity 
fluctuations for strip T5, normalized to the mean value fi and 
standard deviation a, in order to allow comparisons with a 
Gaussian distribution (solid red line). The PDF of the spa- 
tial variations appears to be more peaked compared to the 
Gaussian distribution. A measure of the "peakedness" of a 
probability distribution is given by the kurtosis parameter, k, 
defined as k =< SI >^/<t'', where < SI is the fourth mo- 
ment probability distribution. Since Gaussian distributions 
have k = 3, then the excess kurtosis k is commonly used, which 
is defined as A: = A:- 3. _ 

The excess kurtosis k calculated for the PDF of strip T5 
is 2.44 ± 0.17 indicating a non-Gaussian distribution. Fur- 
ther evidence of non-Gaussian statistics is given by the nor- 
mal probability plot of the cumulative distribution function 
(CDF). This is a quantile-quantile (Q-Q) plot where quanti- 
ties of the observed CDF (y-axis) are plotted against that of 
a normal or Gaussian CDF (x-axis). If the data are normal 
distributed then the normal probability plot of the CDF will 
be linear, while other distribution types will introduce curva- 
ture in the plot. The left panel of Figure |5] shows the CDF of 
strip T5, which does not follow the theoretical function ex- 
pected for a Gaussian distribution (red dot-dashed line). The 
temporal fluctuations are of different character: they are more 
closely described by a Gaussian distribution function. The 
right panel of Figure |4] shows the PDF of the intensity fluc- 
tuations for square D which has k = 0.78 ±0.16, indicating 
statistics very close to Gaussian. Furthermore, the Q-Q plot 



in the right panel of Figure |5] confirms a distribution nearly 
Gaussian for the temporal fluctuations. The quasi-normal dis- 
tribution of the fluctuations in the time domain may again be 
a consequence of the cadence of the observations as discussed 
in the previous section. 

5. STATISTICAL SCALING PROPERTIES OF FINITE 
RANGE TURBULENCE 

A key property of turbulence is that it can be characterized 
and quantified in a robust and reproducible way in terms of 
the ensemble averaged statistical properties of fluctuations. 
We can access to the statistical scaling properties of a spatial 
series, fix), along a given direction x, by constructing differ- 
ences S f with increment L: 

Sf(x,L) = f(x+L)-f(x) (1) 

on the spatial field. Generalized structure functions (GSFs) 
are a powerful tool to test for statistical scaling and are defined 

as: 

/oo 
\Sf\"P(Sf,L)diSf) , (2) 
oo 

where the angular brackets indicate an ensemble average over 
X, implying an assumption of approximate statistical homo- 
geneity. Fully developed inertial range turbulence in an in- 
finite medium exhibits the following scaling for the pth mo- 
ment of the GSF: 

S''(L)^L'^^"\ (3) 

where the ((p) are the scaUng exponents, which are generally 
a non linear function of p. 

For the special case of statistical self-similar (fractal) pro- 
cesses one finds a linear form of ((p) in p, such that: 

ap)=pH, (4) 

where H is the Hurst exponent. 

In fluid turbule nce, we antic ipate intermittency, that is ((p) 
is quadratic in p ( lFrischlll995h . Determining the precise ((p) 
is central to testing turbulence theories. Since we do not have 
measurements in-situ here, we cannot directly compare our 
observed ((p) value with predictions of turbulence theories. 
However we can test whether the ({p) that we observe are 
non-linear with p, consistent with a multifractal, intermittent 
flow and we discuss this in the next section. 

First we will focus upon the direct observations of fluctu- 
ations in the spatial field as these capture non-trivial correla- 
tions in the fluctuations in the flow. The left panels of Figure 
|6]show log-log plots of the averaged 3rd moment of the GSF, 
< ^3 >, versus L,iam (top panel) for strips Tl to T5 and versus 
Liong (bottom panel) for strips LI to L6, where L,rans and Lio„g 
identify the pixel increments of the fluctuations SI(L) along 
the transverse and longitudinal direction respectively. Recall 
that < ^3 > refers to the average over the structure functions 
calculated for each of the 10 rows (or columns) forming a sin- 
gle strip. The structure function analysis provides a measure- 
ment of the correlation of the fluctuations with length scale L. 
The increase of the GSFs with L in the left plots of Figure |6] 
thus suggests that the spatial intensity fluctuations of the QP 
are highly correlated. This is a signature of the presence of co- 
herent structures in the flow. In particular, the intensity fluctu- 
ations in the longitudinal direction (bottom-left panel) reveal a 
correlation over a broader range of spatial scales as the coher- 
ent structure detected are associated to the up and down flows 
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of the QP, which move along the vertical direction; this is the 
longitudinal direction in whic h we expect to see the strongest 
correlation in a turbulent flow (Frisch!'1995'). Along the trans- 
verse direction (top panel of Figure |6l), the curves exhibit a 
knee within a range of length scales of 0.9-2 Mm (dashed 
black lines). These "break points" delimit the crossover be- 
tween the small-scale turbulence and the large-scale coherent 
structures; the above length scales have been attributed to the 
typical distances between the dark up-flows and seem to be 
in good agreement with the multi-mode regime of Rayleigh- 
Taylor instability (Rvutova et al. 2010). 

6. EVIDENCE FOR ESS AND MULTIFRACTAL 
SCALING 

The GSFs shown in the left panels of Figure |6] clearly do 
not follow the power-law scaling of Equation (O. Corrections 
to this equation indeed have to be taken into account for real 
turbulent fl ows for w hich finite range turbulence effects may 
arise ( Dub rullell200Q : Bershadskii 2007); either when turbu- 
lence is not completely evolved (low Reynolds number), the 
data sets are of finite size (realistic cases) or the system is 
bounded, then symmetries in the flow are broken, and the sim- 
ilarity is lost. Nevertheless, a generalized similarity or ESS 
has been observed, which suggests a generalized scaling for 
the pth moment of the GSF by replacing L in Equation Q by 
an initially unknown function G(L), such that 

5^(L) - GiL)'^'-''^ . (5) 

This arises directly from ESS jBenzi et al.l 1 1 9931; 
iCarbone et al.l Q996). Comparing structure functions of 
different order p and q, we can then write: 

5^(L)= [m]^*"'^^^" . (6) 

A log-log plot of S'' versus S'' will therefore give the ratio of 
the respective scaling exponents, Cip)/CiQ)- 

The middle panels in Figure|6]show the ESS in logarithmic 
scale of the GSF for p=2 and q=3. These are straight lines 
on the log-log of Figure |6] thus confirming that Equation (|6]l 
holds. The gradient of such plots in the inertial range pro- 
vides a measurement of the ratio C(2)/C(3). Departures of the 
curves from a linear behaviour occur for length scales outside 
the inertial range and are associated with large-scale coher- 
ent structures in the flow. Finally, the right panels in Figure 
|6]show C(2)/C(3) for strips Tl to T5 along the transverse di- 
rection (top panel) and for strips LI to L5 in the longitudinal 
direction (bottom panel). The error bars provide an estimate 
of the uncertainty in the gradients of the fitted lines in the 
inertial range. The ratio C(2)/C(3) appears to be roughly con- 
stant across all the strips and, more interestingly, differs from 
the value that one would expect if ((p) was linear in p, i.e. 
C(2)/C(3) = 2H/3H - 0.66 (see Equation ©). The ratios of 
the scaling exponents found for all the strips are therefore con- 
sistent with a non-linear form of the scaling exponent C(p)- 
This is a signature of the multifractal nature of this system 
which indicates intermittency within the QP flow. 

The generalized similarity has been tested exphcitly in the 
inertia l range of solar wind turbulence by e.g. IChapman et al.l 
(120091) who formalized Equation (|5]l as follows: 

SP(L) = [SP(Lo)]G(L/Lo)'^^P'> , (7) 

where Lq is some characteristic length-scale of the flow. We 
finally test this generahzed scaling for the intensity fluctua- 



tions of strip T5. In Figure |7] we plot the 3rd moment of 
the structure function, S3, normalized to a value Lq ^ 0.54 
Mm against L/Lo (in logarithmic axes) for 7 consecutive time 
intervals separated by AT = 1.12 min and starting at fo= 
01:10:31 UT. The choice of the value for the parameter Lq 
arises from the characteristic width (on average) of the up- 
flows. The collapse of all the GSFs onto each other within 
the inertial range indicates the existence of a single scaling 
function G{L/Lq). The overlapping of the various GSFs in 
Figure |7]breaks where the effects of large-scale structures be- 
come important. This break point for the strip T5 occurs at a 
length scale of ^^2.8 Mm, which corresponds to the width of 
the large bright structure shown in Figure[T]on the right of the 
prominence and crossed by strip T5 and the length scale of 
transition to coherent structures in Figure |6] (top-left panel). 

7. CONCLUSIONS 

We performed the first qualitative test for in-situ turbulence 
in a QP observed by Hinode/SOT. We analysed the statistical 
properties of the spatio-temporal intensity fluctuations associ- 
ated with the imaged QP from the prospective of a finite sized 
turbulent system. We found the following: 

In space: 

1 . The PSDs of the intensity measurements in the space do- 
main exhibit power law scaling suggestive of non-trivial dy- 
namics. 

2. The PDFs of the intensity fluctuations are described by 
non-Gaussian statistics consistent with small-scale MHD tur- 
bulence. 

3. The GSFs of the intensity fluctuations suggest a gener- 
alized scaling for the structure functions with a dependence 
on a function G{L). They also reveal a high degree of cor- 
relation especially along the longitudinal direction to the bulk 
(driven) flow. Characteristic length scales in the transverse di- 
rection have been detected and associated to the characteristic 
distances between the up-flows. 

4. ESS holds for all the strips considered and it is consistent 
for each direction transverse and longitudinal to the flow as 
a signature of the generalized similarity expected for finite 
range turbulent systems. 

5. The ratio of the scaling exponents C(2)/C(3) is roughly 
constant for all the strips along each direction and its value 
is distinct from 0.66, that is the value expected for a frac- 
tal system. The prominence flow is therefore multifractal in 
character, again consistent with in-situ turbulence. 

6. The intensity fluctuations in the space domain satisfy the 
generalized scaling anticipated by ESS and a scaling function 
G(L/Lq) is observed for different successive time intervals. 

In time: 

7. The PSDs show ~ 1// scaling, consistent with uncorre- 
lated pulses moving past the line of sight of the observations; 
The intensity fluctuations are close to Gaussian distributed. 

The principal aim of this paper has been to explore, for the 
first time, the possibility of discerning the quantitative signa- 
tures of turbulence, namely multifractal or intermittent statis- 
tical scaling, within the flows of a long-lived QP. We have 
shown how tests for non-Gaussianity, multifractality, scaling 
and ESS can be applied in order to fully identify and quan- 
tify statistical properties of turbulent fluctuations. For these 
specific intensity measurements we are restricted to a qualita- 
tive characterization of the fluctuations since the observations 
are integrated along the line of sight rather than in-situ in the 
flow. Despite this constraint, the statistical methods used are 
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Table 1 

Spectral indices 



Domain 


Data 


a ± Aa 




T4 


3.17 ±0.15 


Wave number 


T5 


2.93 ±0.19 




LI 


2.73 ± 0.29 




L2 


2.74 ± 0.37 




A 


1.21 ±0.04 




B 


1.17 ±0.04 


Frequency 


C 


1.29 ±0.04 




D 


1.27 ±0.04 




E 


1.20 ±0.04 



powerful tools to test the hypothesis that in-situ flows are tur- 
bulent. Their application indeed revealed that the statistical 
properties of the intensity fluctuations associated with the QP 
of interest are consistent with a MHD turbulent flow for sys- 
tems of finite size. This is a clear evidence of in-situ evolving 
small-scale turbulence within the prominence flow. 

Since many QPs in the Hinode/SOT database exhibit simi- 
lar dynamics, then this opens up the possibility of using these 
QPs as a 'laboratory for turbulence', to investigate for exam- 
ple finite sized effects on the turbulent flow. The question that 
immediately arises is whether the flow in these prominences 
is more generally found to be turbulent. It would be intriguing 
to determine if or how the presence of turbulence in QPs cor- 
relates with their physical properties. Importantly, turbulence 
is a mechanism by which directed flow is transformed into 
heat. Heating at the loop foo tpoints is known t o drive conden- 
sations at the loop tops (e.g.. lKarpen et alj|2001ft . Rather than 
heating driven by a coronal or a chromospheric reconnection 
process, the evidence of turbulence presented here suggests a 
continuous heating supply that could account for the continu- 
ous formation process of QPs. 
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UK Engineering and Physical Sciences Research Council and 
the Science and Technology Facilities Council on the CFSA 
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Figure 1. QP observed by SOT in the Ca II line on November 30th, 2006 at 01:10:31 UT. Image resolution: 1 pixel ~ 77.22 km on the solar surface. The image 
has been rotated to the horizontal position with respect to the solar Hmb. Intensity levels increase from blue to white. The white grid and 5 squares are shown as 
reference for the analysis in the space domain (transverse strips Tl to T5 and longitudinal strips LI to L5) and in the time domain (squares A, B, C, D and E). 




Figure 2. Top panels: Intensity series I{r) in the space domain for strip T5 (left) and the con'esponding first differences 51 = I(r+L)—I(r) with L=l pixel (right). 
Bottom panels: Intensity series /(;) in the time domain for square D (left) and the con'esponding first differences 51 = I(t + T)—I(t) with t = At = 16.8 sec (right). 
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Figure 3. Log-log plots of the intensity spectra for strips T4, T5, LI and L2 (left) and for squares A, B, C, D, and E (right). All the spectra are shifted in the 
y-direction for clarity. PSDs in the wave number domain (left) reveal two regions with different scaling exponents: kr- (solid line) and kr"" (dashed line), while 
in the frequency domain (right) the PSDs show a single scaling, with spectral index a (dashed line). All the a values are given in Table[T] 




(5 I - <5 l>)/ a (5 I - <5 l>)/ a 



Figure 4. PDFs of the intensity fluctuations in space SI = I(r+L)-I{r) for strip T5 with L = 15 pixel '~ 1 . 16 Mm (left) and in time SI = I{t + T)—I{t) for square 
D with T = 1.12 min (right). Both PDFs are normalized to the mean <5/> and the variance a of the intensity. Red solid lines are Gaussian PDFs with ^ = () and 
o-=l. 
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Figure 5. Normal probability plots (Q-Q plots against a Gaussian) of the CDFs of the intensity fluctuations for strip T5 in space (left) and square D in time 
(right). Length and time scales are those in Figure|4] Dashed red lines refer to the probability expected for a Gaussian distribution. 
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Figure 6. Left panels: Log-log plot of < > versus Ltmns for the intensity fluctuations along the transverse direction for strips Tl to T5 (top). Dashed lines 
delimit the range of scales where the curve exhibit a knee. The bottom panel shows the log-log plot of < 53 > vs. L/ong for strips LI to L5 along the longitudinal 
direction to the flow. All the curves are shifted in the y-direction for clarity. Middle panels: Log-log plots of S2 against for all the strips considered showing 
evidence of ESS. Dashed lines correspond to the linear regression fits across the range Lnrms = 0.2-5 Mm for the transverse strips (top) and Lio„g = 0.2-10 Mm 
for strips LI to L5 (bottom). All the curves are shifted in the y-direction for clarity. Right panels: Gradients of the linear fits shown in the respective middle 
panels, which correspond to the ratio of the scaling exponents ^(2)/c^(3). Notice that the ratios are roughly constant for each direction and different from the 
value 0.66 for both transverse (top) and longitudinal (bottom) strips. Values of <^(2)/c^(3) ^0.66 indicate that C,(p) is non-linear in p and therefore the prominence 
flow shows a multifractal character. 
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Figure 7. Log-log plot of < 53(Lfra„j) > / < SiiLo) > versus Ltrans/LQ with Lo= 0.54 Mm for strip T5 for seven different frames of the dataset separated by AT 
= 1.12 min. fo is the time of the observation corresponding to the first frame of the dataset. Notice that the various structure functions superimpose up to scale L 
~ 2.8 Mm. 



